The inverse spectral problem for first order systems 

on the half line 

Matthias Lesch and Mark Malamud 

Dedicated to the memory of M. G. Krein on the occasion of the 90th anniversary of his birth 

Abstract. On the half hne [0, cxd) we study first order differential operators 
of the form 

where B := 'jj ^ , -Bi , -B2 £ M(n, C) are self-adjoint positive def- 

inite matrices and Q : R+ —> M(2n,C), R+ := [0, 00), is a continuous 
self-adjoint off-diagonal matrix function. 

We determine the self-adjoint boundary conditions for these operators. 
We prove that for each such boundary value problem there exists a unique 
matrix spectral function a and a generalized Fourier transform which diag- 
onalizes the corresponding operator in L^(R,C). 

We give necessary and sufficient conditions for a matrix function a to 
be the spectral measure of a matrix potential Q. Moreover we present a 
procedure based on a Gelfand-Levitan type equation for the determination 
of Q from a. Our results generalize earlier results of M. Gasymov and B. 
Levitan. 

We apply our results to show the existence of 2n x 2n Dirac systems with 
purely absolute continuous, purely singular continuous and purely discrete 
spectrum of multiplicity p, where 1 < p < n is arbitrary. 
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1. Introduction 

We consider the differential operator 



where 

'Bi 



B :-- 



-B2 



Bi,B2 € M(n, C) are self-adjoint positive definite matrices and Q : — > 
M(2n,C), M+ := [0, oo), is a continuous self-adjoint matrix function. If Bi = 
B2 = In then (|1.1D is a Dirac operator. 

It turns out that the operator ( |1.1[ ) subject to the boundary condition 

/2(0) = i//i(0) with Bi=H*B2H (1.2) 

generates a self-adjoint extension Lh of the minimal operator corresponding 
to L. Here, /i(0),/2(0) denote the first resp. last n components of the vector 
/(O). 

Let Y(x, A) be the 2n x n matrix solution of the initial value problem 

LY = XY, y(0,A) = Q). (1.3) 

We will prove that there exists a unique increasing right -continuous n x n 
matrix function C7(A),A G M, (spectral function or spectral measure) such that 
we have the symbolic identity 



Y{x, X)da{X)Y{t, Xy = 5{x - i)l2„. (1.4) 

The main purpose of this paper is to investigate the inverse spectral problem 
for the operator Lh- This means to find necessary and sufficient conditions for 
a n X n matrix function a to be the spectral function of the boundary value 



problem (IT), (|TT 



For a Sturm-Liouville operator this problem has been posed and completely 



solved by I. Gelfand and B. Levitan in the well-known paper |10| (see also |16 



23|, |26|). Later on M. Gasymov and B. Levitan proved similar results for 2x2 



Dirac systems Q, l23|^hap. 12] (see also g and 0). 

We note that in |23|, Chap. 12] the determination of a potential Q with 
prescribed spectral function a is incomplete. The self-adjointness of Q is not 
proved. 
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The paper is organized as follows. In Section 2 we present some auxiliary 
results. In particular we prove the self-adjointness of the operator Lh- 
In Section 3 we introduce the generalized Fourier transform 



oo 



(^/f,Q/)(A) := / Y{x,\rf{x)dx 







(see ( ^.41) ) for / G L^Q„^p(M+, C^") and establish the existence of an n x n matrix 
(spectral) measure a such that the Parseval equality 

if, 9) C^r.) = {^H,Qf, ^H,Qg)Ll(R), (0') 

which is equivalent to ( |1.4| ), holds. In the proof we follow Krein's method 
of directing functionals ||14| |, |1T5| ]. Moreover we show that ^h,q is a unitary 
transformation from L^(M-|-,C^") onto L'^i^) which diagonalizes the operator 
Lh- Namely, ,^h,qLh-^h]q = A where A : L^(M) Ll{W) denotes the 
multiplication operator by the function A ^ A. Similar results (with similar 
proofs) hold for Sturm-Liouville operators as well as for higher order differential 
operators. 

In Section 4 we introduce (under the additional assumptions on B) a tri- 
angular transformation operator I + K and present a sketch of proof of the 
representation Y{-,X) = {{I + K)eo){- , X) where 60(3;, A) is the solution of (0) 
with (3 = 0. Then we derive the linear Gelfand-Levitan equation 

F{x,t) + K{x,t) + K{x,s)F{s,t)ds = 0, x > t, (1.5) 
^0 

with F{x,t) defined by ( 4.34| ). F is the analog of the so-called transition func- 



tion (cf. [[16|| ). We present two proofs of (1^). The proof after Theorem 4^ 



is close to the proofs in |10| and [23, Chap. 12]. The second one is rela- 



tively short. It is based on simple identities for kernels of Volterra operators 
(see ( 4.17| )- ([4.23| )). In Proposition 4^ we derive two representations ( 4.31 ) and 



( [4. 34]) for F{x,t) which easily imply (|1.5| ). In other words, this proof derives 
the linear equation (^]^) directly from the nonlinear Gelfand-Levitan equation 
(1^. This proof seems to be new and is essential in the sequel. 



Furthermore, in Section 5 we solve the inverse problem (Theorem ^.2| ). 
Namely, starting with the transition matrix function F[x, t) of the form (^.1|') we 
prove the existence of the unique solution K{x,t) of (^]^). Conversely, starting 
with K{x, t) we determine the matrix potential Q{x) = iBK{x, x) — iK{x, x)B 
and we prove that y(-. A) := ((/ -|- -ftr)eo)(-, A) satisfies the initial value problem 
(P). 



We present several criteria for the prerequisites of Theorem 5.2 to hold. 

Finally, in Section 6 we present some generalizations and improvements of 
the main result. The degenerate Gelfand-Levitan equation is also considered 
here. We point out that we have obtained a sufficient condition for an increasing 
matrix function a to be the spectral function of the operator L. In the special 
case that B = (Ail„, — A2-^«) (or more generally for the class (Tb), cf. Section 
4 ) our conditions are also necessary. Finally, we prove the existence of 2n x 2n 
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systems with purely absolute continuous, purely singular continuous, and purely 
discrete spectrum of any given multiplicity p, 1 < p < n. 

In conclusion we mention some recent publications close to our work. D. 
Alpay and I. Gohberg have constructed some explicit formulas for the 

matrix potential of a Dirac system ( [L.l| ) from the rational spectral function. 
Their approach is based on the results of minimal factorizations and realizations 
of matrix functions @. 

A new approach to inverse spectral problems for one-dimensional 
Schrodinger operators with partial information on the potential as well as 
to different kinds of uniqueness problems on the half-line has been recently 
proposed by F. Gesztesy and B. Simon (see p2| , [p!3|] and references therein). 



Furthermore, we mention the recent paper F. Gesztesy and H. Holden [11] on 
trace formulas for Schrodinger-type operators. 

The results of this paper have been announced in [20], a preliminary version 
of this paper has been published in ]llt 
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2. Preliminaries 

We consider again the operator ( [I.l| ) from the introduction. In the sequel 
for a vector v G C^" the vectors vi,V2 G C" will denote the first resp. last n 
components of v. In this paper scalar products will be antilinear in the first 
and linear in the second argument. This is necessary since we will be dealing 
with vector measures (see ( |3.1D below). 

L is a formally self-adjoint operator acting on i/^Qj^p((0, oo), C^") C 
L2(M+,C2"). We denote by L* the adjoint of L in , C"^"") ■ To obtain 

self-adjoint extensions we impose boundary conditions of the form 

H2f2{0) = HihiO). (2.1) 

Here, Hi,H2 £ M(n,C) and /i(0),/2(0) G C" denote the first n resp. last n 
components of /(O), where / G iJcVp(R+, C^"). 

Proposition 2.1. Let Lhi,H2 be the operator L* restricted to the domain 

&{Lh„h,) ■■= {/ e 1 i^2/2(0) = i^i/i(o)}. 

Then the operator Lh^.H2 ^•s self-adjoint iff the matrices Hi,H2 are invertible 
and Bi = H*B2H, where H := H:^^Hi. 
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Consequently we have Lhi,H2 — ^h^^Hi i ~ ^h.i- From now on we will 
denote Lh,! by Lh and we will write the boundary condition always in the form 

f2{0) = Hfi{0). (2.2) 

Proof. Since Q is continuous we have &{L*) C i7iQ^(M+, C^"). Now choose 
a sequence of functions Xm G Co°(I^) with the following properties: 

(i) Xm\{-oo,m] = 1, 

(ii) < < 1, 
(hi) \x'm\ < 7Pl- 

If / G ^{L*) then Xmf ^ / in L2(R+,C2") and 

Lxmf = B-^xU + Xn^Lf ^ Lf (2.3) 

in L2(R+,C2"). Thus Xmf ^ / in ^{L*). 
For f,ge &{L*) we then find 

{L*f,g) - {f,L*g) = \\m.k^^\\m.i^^{L*Xkf,Xia) - {Xkf,L*xig) 

= -i limfc^oo linii^oo < Bxkf{0),Xl9{0) >c2n= -i < Bf{0),g{0) >c2n . 

(2.4) 

Hence, g G ^{L^j^^h^) ah / e ^(L^^./fJ 

=< 5/(0), 5(0) >c2n . (2.5) 

This shows that any self-adjoint extension of L is given by a Lagrangian 
subspace V of the symplectic vector space C^" with symplectic form 

u){v,'w) :=< Bv,w > = < i?it7i,tt;i > — < B2V2,W2 > ■ 

Lagrangian means that diml/ = n and u!\V = 0. The domain of such an 
extension then is 

{/e^(L*)|/(0)en- 

Now let y be a Lagrangian subspace of C^" = ® C". We denote by 
7ri,7r2 the orthogonal projections onto the first resp. second factor. Since the 
symplectic form uj is positive resp. negative definite on kervri resp. ker7r2 and 
since dim V = n the maps tti , tt2 restricted to V are isomorphisms 

^i-.V^Cl, 1t2--V^C1. (2.6) 

Hence V = {(x,^2 ° t^i^x) \ x G CI). Put H := ^2 o \ Then uj\V = 
immediately implies Bi = H*B2H. This proves the proposition. □ 

Remark 2.2. 1. The previous proposition shows that the deficiency indices 
n±{L) are equal to n, i.e. n±{L) = n. This means that at infinity we do not 
have to impose a boundary condition. Thus infinity is always in the 'limit point 
case', which essentially distinguishes first order systems from Sturm-Liouville 
operators and higher order differential operators ([^, |^]). 

2. For scalar Dirac systems (n = Bi = B2 = 1) another proof of Proposition 



2.1 has been obtained earlier by B.M. Levitan p3. Theorem 8.6.1]. 
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The present proof is adapted from the standard proof of the essential self- 



adjointness of Dirac operators on complete manifolds (see e.g. |18, Theorem 
II.5.7]). 



3. At the same time as our preprint |19| the paper Sakhnovich |30| ap- 
peared. Following Levitan's method he obtained some sufficient conditions for 
a canonical system to be selfadjoint. This is a system 

J'-^ = ^Hi.)yi.,X), J=[\ ^-), (2.7) 
where H{x) is a continuous nonnegative 2n x 2n matrix function. The method 



of proof of Proposition 2.1 can be extended to arbitrary first order systems, 



in particular to generalize the recent result from |30| for canonical systems. 
Details will be given in a subsequent publication. 

4. Another proof of the previous proposition could be given using the 
uniqueness of the solution of the Goursat problem for the hyperbolic system 
^ = ±iL*jjU in M^. This method (see |^) was also used to prove the essential 
self-adjointness of all powers of the Dirac operator on a complete manifold (cf. 
Q). Sakhnovich's result [30] mentioned before also follows from the hyperbolic 



system method. 

For the problem considered here we prefered to present an elementary direct 
proof. 

From now on we will assume 

Bi = H*B2H. (2.8) 

Note that this implies that H is invertible. 

We first discuss in some detail the case (5 = 0. Let A E M(n, C) be a 
positive definite matrix. Then we put for / E L^(M, C") 



^a/(A) := / e-*^" ^''f{x)dx. (2.9) 
Then we have for f,g£ L^(IR, C") the Parseval equality 

if, 9) = ^ [ {^AmrA-\^A9){\)d\. (2.10) 

To prove ( p. 10 ) we may assume A to be diagonal, i.e. A = diag(ai, . . . ,a„), 
because iiA = UAU* with a unitary matrix U then (^a/)(A) = U{,^^U* f){\). 
Now ^AfW = {'"^ fji^l CLj))j=i,... ,n and ( 2.10| ) follows easily from the Parseval 



equality for the Fourier transform. 
Now let 



eo(.,A):=e--^^fiUf (2.11) 



and put 

^H,of{\):= eo{x,Xyf{x)dx = ^Bji{X)+H*^Bj2{-X), (2.12) 
Jo 

where fj denotes the extension by of fj to M. 
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li f,ge L2(M+,C2") then the integrals 

^ (2.13) 

JR 

are sums of scalar products of the form 



/ 

JR 



^^(-A)^V(A)dA, (2.14) 



where 9?, "0 ^ L^(R-(_). These scalar products vanish and hence we end up with 
the Parseval equality in the case of Q = 



2tt 



+ ^ [ ^Bj2{-XrHB^^H*^BM-X)dX 

= {f,g), (2.15) 

in view of Q and (l2lc| ). 

3. The spectral measure 

In this section we prove the existence of a spectral measure function for the 



self-adjoint operator Lh based on Krein's method of directing functionals |14 



151. For the convenience of the reader we recall Krein's result. 



Definition 3.1 (]14|, p!5|]). Let A be a symmetric operator in a separable 
Hilbert space H and let S be a dense linear subspace of H containing ^{A). 

The system of linear functionals defined on E and depending on 

A G M is called a directing system of functionals for yl in if the following 
three conditions are fulfilled: 

1. ^ j{f; X), j = 1, ...,p, is an analytic function of A G M, for each f ^ E; 

2. the functionals $j(-; Aq) are linearly independent for some Aq G M; 

3. for each /qGE and Aq G M the equation Ag — X^g = /o has a solution 
in E if and only if 

^i(/o;Ao) = foralll<j<p. 



Theorem 3.2 ( |14f| , |15i ). Let A be a symmetric operator in H with 
S'{A) C E C H which has a directing system of functionals {^j{-;X)}^ in 
E. Then 

1. there exists an increasing p x p matrix function a {X) = icjkiX))^ j^^^ such 
that the equality 

P r 

i9J)= E / ^ji9;X)Mf;X)dajk{X) 
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holds for each f,g £ E. 

2. // a is normalized by requiring it to be right-continuous with cr(0) = 
then it is unique if and only ifn+{A) = n^{A), where n±{A) := dimker(yl* =pz) 
denote the deficiency indices of A. 

Definition 3.3. Let C7(A) = {o'ij{X))fj^i be an increasing n x n matrix 
function. On the space Co(M, C") of continuous C^-valued functions with com- 
pact support we introduce the scalar product 

if,9)Ll ■■= [ f{\rda{\)g{X) := V / W}g,{X)dai,{X). (3.1) 
We denote by L'^{^) (cf. p7|]) the Hilbert space completion of this space. 



Remark 3.4. From now on we will consider - without saying this explicitly 
- only right-continuous n x n matrix functions which map to the 0-matrix. 
Such a function a is determined by its corresponding matrix measure da. 

We turn to general Q. For future reference we state the boundary value 
problem for L: 

Lf = Xf, /2(0) =F/i(0), where Bi = H*B2H. (3.2) 

Proposition 3.5. Let Y : M+ x C ^ M(2n x n,C) be the unique solution 
of the initial value problem 

LY{x, X) = XY{x, A), y(0, A) = ^) . (3.3) 

Then: 

1. There exists an increasing nxn matrix function C7(A),A G M, (spectral 
function) such that the map 



poo 

■^H,Q : L2„^p(M+,C2") 9 / ^ (^h,q/)(A) := F(A) := / Y[x,Xrf{x)dx 

Jo 

(3.4) 

extends by continuity to an isometric transformation from L^(M_|_, C^") into the 
space L^(M), i.e. for f,g e L^(M+,C^"') we have the Parseval equation 

f*{t)g{t)dt = [ F*iX)da{X)G{X) (3.5) 

JR 

with F,G being the ^h,q "transforms of f,g. 

2. If a is normalized by requiring it to be right-continuous with (t(0) = 
then it is unique. 

Proof. 1. Let b G M+ be a fixed point and let Lb be the operator L* 
restricted to the domain 

^(Lb) = {/ G H\[0,b],C^n I /2(0) = if/i(0),/i(6) = /2(6) = 0}. (3.6) 
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It is clear that Lf, is a symmetric operator and = (Lb)* is a restriction of L* 
to the domain 

^(L^) = {/ G H\[0,b],C^n I /2(0) = ^/i(0)}. (3.7) 

We consider ^(L^) as a subset of (M.^ , C"^"^) identifying each function 
/ G Si {Lb) with its continuation by zero to M^.. 

Since is a regular differential operator on a finite interval, each A G C 
is a regular type point for L^-, i.e. \\{Lb — A)/|| > e||/|| for all / G Si{Lb) with 
some e > 0. In particular, — X has closed range. Hence, for a fixed A G M 
and / G L2([0, 6]; C^") the equation 

Lg-Xg = f, A G M, (3.8) 

has a solution 5 G L^([0, C^") if and only if / is orthogonal to the kernel 
ker(L* - A), that is if JoY*{x, X)f{x)dx = 0. 

Denoting by Yi the i-th column of Y ^ on rewrites the last equation as 



$,(/;A) := / <Y,{x,X)J{x)>dx 
Jo 

= / Yji{x, X)fj{x)dx = 0, l<i<n 



i=i 

It is clear that the functionals on L^Qjjjp(M+, C^"), defined by the left-hand 
side of (|3.9D, are linearly independent and holomorphic in A G M. Thus the 



conditions 1. and 2. of Definition 3.1 are satisfied. Since E := L^^j^^ (M-|_,C^") 



is dense in L (M-|_,C ) the functionals $j(/, A) thus form a directing sys- 
tem of functionals for the operator A := Lh \ &{Lh) H L1^^^{E.j^,£?'^). 
By Krein's Theorem |3.2| there exists ct(A) such that (|3.5| ) holds for arbitrary 

-^comp(I^+) C "). 



2. In view of Proposition 2.1 the operator A = Lh \ &{Lh) n 



L^ojjjp(M+, C^") is essentially selfadjoint and consequently n+{A) = n-{A) = 0. 
Thus the uniqueness of <t(A) follows from the assertion 2. of Krein's theorem 



3.2. □ 



Remark 3.6. 1. The Parseval identity may be symbolically rewritten as 
Y{x, X)da{X)Y{t, Xy = 6{x - t)hn. (P') 



To obtain (U') from (|j) it suffices to set in (^) /(^ = 4(^) ® 
^ii — ^tiO ^ hj ^ 2n and to note that (^hq/)(A) = 

Yix, Xfei, {^H,Q9){X) = Y{t, A)*e,. 



2. Another proof of Proposition 3.5 based on the approximation method 



proposed independently by B.M. Levitan [^^, Chap. 8] and N. Levinson 0, 
Chap. 9] was given in the preliminary version of this paper [pT " 



For convenience we denote the extension of '^h,q to L^(M-|_,C^") by the 
same letter. Next we prove the surjectivity of '^h,q- 
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Theorem 3.7. Under the assumptions of Proposition |3.5| the mapping 
^H,Q is surjective, that is ^h,q raaps (M^- , C^" ) onto L^(M). 

Proof. So far we have proved that ^h,q ■ L2(M+,C2") — > ^^(R) is an 
isometry. To prove surjectivity we mimick the proof of [|^, Sec. 9.3] for second 
order operators. 

Note first that for / e Si[Lh) we have 

{^H,QLHf){\) = -Y{0, Xy^BfiO) + A(^h,q/)(A) = X{^H,Qf){X) (3.10) 



since in view of (^) and (|3J) Y{0, A)*5/(0) = 0. 

For / G Ll^^^(R+, C'^'')n&{LH) formula (|3T0| ) follows from integration by 
parts. For arbitrary / G &{Lh) it follows from the fact that Ll^^^i^j^,'^?^) n 
2i{Lh) is a core for Lh- The latter follows from the proof of Propo 

Next we construct the adjoint of '^h,q- we put for g G -^^^,0011 



^Hg){x) := / Yix,\)daiX)giX). (3.11) 



Then for /GL2„^p(M+,C2") 

^Hg,f)L2{R+,c2'-) = I {^Hg)ix)*fix)dx 

g{X)*da{X)Y{x,X)*f{x)dx (3.12) 
From the estimate 

K^HS', /)l2{R+,C2")| = \ {9^'^H,Qf)Ll{B.)\ < IIS'IIl2(ir;)||/||l2(]k_^^c2-) (3.13) 

we infer that '^h extends by continuity for L^(IR). Moreover, it equals the 
adjoint of ^h,q, i-e. 

Since ^h,q is an isometry it remains to prove injectivity of q- 
It follows from ( 3.10| ) and Proposition |2.l| 

^h,q{Lh-0-^ = {^-C)-^^h,q (3.15) 

ioT C = y + ie Here A : Ll{W) L2(M), {Kg){X) := Xg{X) denotes_the 

operator of multiplication by A. Therefore {Lh — C)~^'^hq = '^hq{^ ~ C)~^ 
and hence we have the implication: 

^> E ker Q ^ (A - C)"^^' € ker g for all C G C \ M. (3.16) 

We put 

Y{x,X) ■= [ Y{t,X)dt. (3.17) 
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Note that the i-th row Yi{x, A) of Y satisfies 

Yi{x, Xy = r Y{t, XYadt = (^//.q(1[o,x] ® ei))(A), (3.18) 

JO 

where Cj denotes the i~th unit vector in 
Now let $ e ker^^ Q. 



-<2n 



In particular Yi{x,-)* G L'^{^) and thus in view of and (|3l6|) we have 

for X > 0, e > 0, e M 



r'^) = ffi* w. ,L , (3-19) 



thus 



0= / ^(^'A)- ^da(A)$(A). (3.20) 

Since l^(a;, •)(icT(-)<I>(-) is L} the dominated convergence theorem implies for 
a,/? G M 

= \im [ Y{x,X)- ^^da{X)^X)du 

= I y(x,A)lim / %: ^di/(icT(A)$(A) (3.21) 

Jr ^-"OJa {X-uy + s'^ 

= vr / Y{x,X)da{X)<^>{X). 



Differentiating by x and putting x = yields for a, /3 G 

r-/3 



= / Y{0,X)da{X)^X). 

J a 



Since 1^(0, A) = (^) and H is invertible we have 

da(A)$(A) = 



for all a, /3 G M. This imphes $ = in Ll{R). 

□ 

Remark 3.8. 1. We note that another proof of the uniqueness of the spec- 



tral function in Proposition 3.5 can be given using Theorem 3.7. To prove 
the uniqueness statement we assume we had another increasing right continu- 
ous n X n matrix function g, q{0) = 0, such that ^h,q is a unitary transfor- 
mation from L^(M+,C^") onto L^(M). Then in view of ( p.5| ) we have for all 
F,GGL2„^p(M,C-) 

F{Xyda{X)G{X) = [ F(A)*d^(A)G(A), 



and hence the two Radon vector measures da and dg coincide. By the right- 
continuity and the normalization g{0) = cr{0) = this implies a = g. 
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2. For n = 1 Proposition p.5| follows from |15, Theorem 4]. We also note 
that a generalization of Krein's theorem to the case n > 1 may be obtained by 
a slight modification of the proof of Proposition 

3. I n |29| , Chap. 3] the existence of the spectral function for a canonical 
system (|2.7|) is stated. For nonsingular Hamiltonians this fact follows from 



Krein's Theorem |3.2| in just the same way as Proposition |3^ 



We note also that for a singular Hamiltonian similar results may be obtained 



by the corresponding generalization of Krein's Theorem for linear relations 



Example 3.9. ( 2.15| ) shows that in the case Q = we can choose for a the 



function (To(A) := 2^B^ A. 

4. Transformation operator and Gelfand Levitan equation 



1. We present a special case of |24, Theorem 7.1], (see also |25, Theorem 



1.2]). In the sequel we assume -B to be a diagonal matrix, which can be achieved 
by conjugating L with an appropriate unitary matrix. 
Let 

B = diag(AiI„,, . . . , A.r/„J, 

r I 1 (4-1) 

ni = minjnj | 1 < z < rj, ni + n2 + ■ ■ ■ + Ur = 2n. 
Furthermore, we put 

n:={{x,t) eM.'^\0<t<x}. (4.2) 



Theorem 4.1. Let B he as in (|]|) and let Q = iQij)lj=i ■ K+ — > 
M(2n, C) continuous, where Qij denotes the block-matrix decomposition with 

T 

respect to the orthogonal decomposition C^" = © C"'. Moreover, we assume 

i=l 

that Q is off-diagonal, i.e. 

Q.. =0, i = l,...,r. (4.3) 
Let Y he the solution of the equation (|3.3| ) satisfying 

y(0, A) = A := col(^i, . . . ,^r), G M(nj X ni,C), rank^j = ni. (4.4) 
Then there exists a continuous function K : Vt — > M(2n, C) such that we have 



Y{x,X) = Yo{x,X)+ I K{x,t)Yoit,X)dt, (4.5) 

where 







Yo{x,X)=e'^^~'''A 



is the solution of the equation (|3.2[) with Q = and satisfying the same initial 
conditions ( [4.4[ ). 

IfQe Ci(M+,M(2n,C)) then K £ (7^(17, M(2n, C)) and it satisfies 

Bd^K{x, t) + dtK{x, t)B + iQ{x)K{x, t) = 0, (4.6a) 
BK{x,x) - K{x,x)B + iQ{x) = (4.6b) 
K{x,Q)BA = Q. (4.6c) 
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If Q ^ C(M+, M(2n, C)) then K is the generalized continuous solution of ( |4.6| ) . 
Conversely, if K is a (generalized) solution o/(4.6) then Y{x,X) defined by 



( [4 .51 ) is the (generalized) solution of the initial value problem ( |3.3| ). 

Sketch of proof, i) Suppose that K G C^(R+, M(2n, C)) and that for- 



mula ( [4. 51) holds. Substituting (4^) into ( 0.3D and integrating by parts one 
obtains 

[BK{x, x) - K{x, x)B + iQix)] Yq{x, A) + K{x, ^)BYo{<d, A) 

r r n (4-7) 

+ / [Bd^K{x,t)+dtK{x,t)B + iQ{x)K{x,t)\YQ{t,\)dt = 0. 
Jo 



Since yo(0, A) = A does not depend on A one concludes from ( [4.7| ) and the 
Riemann-Lebesgue Lemma that (|4.7D is equivalent to (|4.6| ). Thus in this case 
the representation ( [4. 51) is equivalent to the solvability of the problem (|4.6| ). 

ii) Next we prove the existence of a (not unique) solution of the problem 
( 4.6a )~( [l.6b| ). Let R{x,t) be one of them. Using the block-matrix representa- 



tion R{x,t) = {Rij{x,t))l we rewrite the problem ( 4.6aD -( ^.6b ) as 



XidxRijix,t) + XjdtRij{x,t) = -\/^ Qip{x)Rpj{x, t), 1 < i,j < r, 

p=i 

(4.8) 

Rij{x,x) = -V^{Xi - Aj)"^ Qij{x), 1 < i / J < r. 

(4.9) 



It is clear that the system ( [4.8[) is hyperbolic with real characteristics lij : x = 
kijt + c{kij = XjX~^). Thus, m^} = {0<t<x< oo} we have the incomplete 
characteristic Cauchy problem ( [4.8| ), ( [4.9[ ) with (2n)^ — — ... — scalar 
conditions (4.9). Fixing xq £ M-|- and setting 



^min — Uia,x{kij \ kij G (0, 1), ^ ^ i, j ^ r}, A^max — ^miiT 



we consider the triangle Aabc confined by the lines AB : x = t, AC : x — xo = 
kmint, BC : X — xq = /cmax*- We preserve the notation Q{x) for a continuous 
extension to M of the function Q{x) with the same norm. Furthermore, we 
denote by a and b the abscissas of the points A and B respectively. Now we 
impose the following n1 + . . . + n'^ conditions on the characteristic line AC : 

Rjj{x,{x - l)kmm) = 0, for x£[a,xo], j£l,...,r. (4.10) 

Thus, we arrive at the Goursat problem ([4.8|)-( |4T0|) f or the hyperbolic system 
( [4.^ ) in the triangle Aabc- Integrating the system ( |4.8[ ) along the characteristics 



and using ( |4.9| ), ( |4.10 ) one deduces the system of integral equations 
XiRij{x,t) 



Xi — Xj 



'1 / Yl QipiORpjiC, - x)kij + t)dc, 
Hiji^^t) p=i 

(4.11) 



14 MATTHIAS LESCH AND MARK MALAMUD 

where for brevity it is set Aj(Aj — Xj)~^Qij{^ij{x,t)) = for i = j and 

{XjX - \it){\j - Xi)-^, i^j, 



a + {1 — a){x — t), 



For Q G C^{M.,M{n,C) the system (4.11) is equivalent to the Goursat problem 
([4.^)-( O!0| ). The solvability (and uniqueness) of the solution of ( 4.11| ) is proved 
by the method of successive approximations. 

For Q £ C(M, M(n, C)) \ (M, M(n, C) we understand the solution of (gj)- 
( [4.10D as a solution of ( |4.11D . 

iii) To finish the proof, starting with the solution R(x,t) of the Goursat 
problem (|4.8D-(4.10) we introduce a convolution operator 



^{x - t)f{t)dt 



with ^{x) = diag(<I>i(a;), . . . ,^r{x)) being a block-diagonal 2n x 2n matrix 
function, consisting of Uj x Uj blocks and define the operator K by the 
equality I + K = {I + R){I + It is clear that X is a Volterra operator with 
the kernel 

rx 

K{x,t) = R{x,t) + ^{x-t)+ R{x,s)^{s -t)ds. (4.12) 

Jo 

Since the operator I + R intertwines the restrictions Lq and —iB Dq of the 
operators L and -iB ® D onto {/ G H^{[^, 1], C^^^) | /(O) = 0}, that is Lq{I + 
R) = [I + R){—iB (8) -Do)) so is / + K. This fact amounts to saying that K{x, t) 
satisfies the problem (|4.6a]) -( ^.6b ). To satisfy the condition ( [4.6c|) it suffices (in 
view of ( 4.12| )) to choose $(x) as the solution of the equation 

^{x)BA+ [ R{x,s)^{s)BAds = -R{x,0)BA. (4.13) 
Jo 

Since rank j4j = ni,2 < j < r, the Volterra equation (|4.13| ) is of the second kind 
and therefore has the unique solution <I> G C([0, oo), M(2n, C)). Thus K{x,t) is 
the required solution of ( [4.6a| )-( [46c| ). □ 



Corollary 4.2. Under the assumptions of Theorem |4.1| let B = 
(Ai/„,A2/n) (that is r = 2.) 

Then there exists a continuous function K : 0, — > M(2n, C) such that we 
have 

Y{x,X) = eo{x,X) + [ K{x,t)eo{t,X)dt, (4.14) 
Jo 

where eo{x,X) was defined in (2.11). 

IfQ€ Ci(M+, M(2n,C)) then K G 0^(17, M(2n, C)) and it satisfies 

Bd^K{x, t) + dtK{x, t)B + iQ{x)K{x, t) = 0, (4.15a) 
BK{x,x) - K{x,x)B + iQ{x) (4.15b) 

r 



K{x,0)B 



H 



0. 



(4.15c) 
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If Q ^ C(M+, M(2n, C)) then K is the generalized continuous solution of ( |4.6| ) . 

Conversely, if K is a (generalized) solution of ( 4.15| ) then Y{x,X) defined 
by ( |4.14 ) is the (generalized) solution of the initial value problem ( |3.3D . This 
last statement holds even for general B of the form ( |1.1[ ) . 

2. We continue with some general remarks about Volterra operators: 

For any continuous matrix function K : — > M(2n, C) we obtain a 
Volterra operator 

Kf{x) := r K{x,t)f{t)dt (4.16) 

JO 

acting on C(M+,C2") or L2([0, a], C^") for any a > 0. By slight abuse of 
notation we will use the same symbol for the operator and its kernel. The 
set of operators I + K with K being a Volterra operator forms a group. The 
operator 

R := {I + K)-^ - I (4.17) 

is again a Volterra operator with continuous kernel R{x,t),t < x. From the 
equation 

I = {I + R){I + K) = {I + K){I + R) (4.18) 

we deduce 

RK = KR = -R-K. (4.19) 

Put 

F := R + R* + RR*. (4.20) 

The kernel of F obviously is 



F{x,t) = < 



R{x,t)+ I R{x,s)R{t,s)*ds, x > t, 
Jo 

R{t,x)*+[ R{x,s)R{t,s)*ds, x<t. 



(4.21) 



Furthermore, using ( 4.19 ) we conclude 

F + K + KF = R + K + R* + RR* + KR + KR* + KRR* = R* (4.22) 
thus we have the " Gelfand-Levitan equation" 

F + K - R* + KF = 0. (4.23) 

Proposition 4.3. Let K : ^ — > M(2n, C) be continuous and let R : Q — > 
M(2n,C) be the continuous kernel of the Volterra operator {I + K)~^ — I. Then 
the function F : — > M(2n, C) defined by ( 4.21| ) satisfies the "Gelfand- 
Levitan equation" 

rx 

F{x,t) + K{x,t) + K{x,s)F{s,t)ds = 0, x > t, (4.24) 
Jo 

/•X 

F{x,t) - R{t,xy + K{x,s)F{s,t)ds = 0, x<t. (4.25) 
Jo 
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Conversely, if Fi : Q — > M(2n, C) is continuous and satisfies (4.24) then 

Fi = F\n. 

Proof. It only remains to prove the assertion about Fi. The difference 
F{x,t) — Fi{x,t) satisfies the equation 

F{x,t) - Fi{x,t) + K{x,s)[F{s,t) - Fi{s,t)]ds = 0, 0<t<x. 
Jo 

For each fixed t £ [0,x] this is a homogeneous Volterra equation of the second 
kind and consequently has only the trivial solution F{x,t) — Fi{x,t) = 0. □ 



We turn back to the system (|3.2|). 

Definition 4.4. We say that the system ( ^ ) (resp. the operator L) be- 
longs to the class (Tb) if for this system there exists a transformation operator. 

This means that the solution Y{x,X) of the initial value problem ( p. 3D ad- 
mits a representation (4.14) with a continuous function K : ^ M(2?i;C). 
Corollary says that the system is of class (Tb) if i? = (Ail„, \2ln)- 



It follows easily from Proposition 4.6 below that for an operator Lu of class 
(Tb) the transformation operator I + K is unique, i.e. the representation ( p.3| ) 
for y(x, A) is unique. 

If the system is of class (Tb) then we denote by K the unique Volterra 
operator with continuous kernel satisfying ( [4.15 ). As before R denotes the 
Volterra operator defined by R := {I + K)~'^ — /. 

In particular we have in view of ( |4.15| ) 

eo{x,\) = {{I + R)Y{-,\)){x)=Y{x,\)+ f R{x,t)Y{t, X)dt. (4.26) 



Lemma 4.5. Let L be of class (Tb). 

1. Let a be the spectral function of the boundary value problem ( |3.2[) and let 
ffGL2„^p(M+,C2«). Put 

Go(A) := i^H,o9)W = / eoix,\yg{x)dx. (4.27) 

Jo 

Then Gq G LI{R) and if 

[ Go{XTda{X)Go{X) = (4.28) 

then g = 0. 

2. We have ^H,o(^Lmp(IR+, C^")) = ^^,,q(L2„^p(R+, C^")). 
Proof. In view of ( |4.26D we have 

[Y{x,X)* + J Y{t,XyR{x,tydt]g{x)dx 

/oo 
R{t,xyg{t)dt]dx, (4.29) 
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hence Go(A) is also the ^/f^g-transform of the function 

g{x) := ((/ + R*)g)[x) = g{x) + / R{t, xrg{t)dt. (4.30) 



Since g G L^ojnp(M+, C^") we also have g G L^omp(J^+) C^"). This shows 
the inclusion #'//,o(icomp(I^+> C^")) C ^H,Q(icomp(I^+. C^")). The converse 
inclusion is proved analogously using (|4.5| ) instead of (|4.26D . 

In view of the Parseval equality (Proposition 3.5) we find 



/ 



Go{Xyda{X)Go{X) = / g{xyg{x)dx 



which by assumption ( [1.28 ) implies 5 = 0. Since g has compact support ( 4.30| ) 
is a Volterra equation and thus g = 0. □ 



Proposition 4.6. Let a be the spectral function of the boundary value prob- 
lem ( [3.2| ) and let (Tq = ^B^^X be the corresponding spectral function for Q = 0. 
We abbreviate S := a — aQ. 

1. Let L be of class (Tb) and let I + R be the transformation operator of 
the form (4.26). Furthermore, let F be the 2n x 2n matrix function defined by 

R{x, t) + [ R{x, s)R{t, s)*ds, x >t>0, 
Jo 



F{x,t) := < 



R{t, x)* + / R{x, s)R{t, s)*ds, 0<x<t. 
^ Jo 



(4.31) 



Then we have for all f,g^ -^comp( 
Fo(A)*dS(A)Go(A) : 



JO 



f{x)*F{x,t)g{t)dxdt, 



(4.32) 



where Fq^Gq denote the c^^Hfl^t fans forms of f,g. 
2. Again assuming L to be of class {Tb) we put 



eo{x,X) := / eQ{t,X)dt. 



Then the function 



Fix. 



t) := / eo(x,A)dS(A)eo(t,A)* 
Jr 



(4.33) 



(4.34) 



exists and has continuous mixed second derivative which coincides with 
F{x,t), i.e. ^^F{x,t) = F{x,t). 

3. Conversely, given any increasing n x n matrix function a put S := 
(7 — (Tq. If the integral ( 4.341 ) exists and has a continuous mixed second derivative 
Fi{x,t) := ^^F{x,t) then O holds for all f,ge L2„^p(M+,C 



2n\ 



dxdt 

instead of F. 



with Fi 
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Remark 4.7. We emphasize that 3. holds for arbitrary L of the form ( p.2| ) 
not necessarily being of class (Tb). 

We note that the identity ( 4.32[ ) characterizes the spectral function of the 
problem (|3.2|). More precisely, if q is an increasing (normalized) n x n matrix 



function such that (4.32) holds with "Eg := g — ao then g = a. 
Indeed from ( 4.321) we infer 

Fo(A)*da(A)Go(A) = / Fo(A)*(i£>(A)Go(A), ^ ^ 

(4.35) 

Fo ■■= ^Hfif, Go := ^H,og, 



for all f,g^ L^^^^p (M_(_ , C^"). By Theorem and Lemma 4^, 2. this implies 
that ( [4.35D holds for ah Fo,Go G Ll{R), in particular it holds for ah Fo,Go G 
C(]R, C") with compact support. Thus the vector measures da,dg and hence 
the right-continuous functions g, a coincide. 

Proof. 1. In view of ( [4.29|) Fq is the g-transform of 

f{x) = f{x)+ R{t,xrf{t)dt, (pO') 

J X 

thus the Parseval equality (|3.5|) gives 

Fo(A)*(fS(A)Go(A) = / Fo(A)*d<T(A)Go(A)-(/,g) = (7,5) -(/,<?) 
(/(x)+/ R{t,xyf{t)dtng{x)+ R{t,xyg{t)dt)dx-{f,g) 



f{x)*F{x,t)g{t)dxdt 

'0 ^0 

by a straightforward calculation. 

2. For a;,t > and /o,5o G C^" we apply 1. with f{u) := l[Q^x]{u)h, g{v) := 
Mo,t]iv)go and find 

/o* / eo{x,X)dE{X)eo{tAr9o = fo T ^ F{u,v)dudvgo, (4.36) 

Ji? JO JO 

which implies the first assertion. 

3. To prove the converse statement we note that now we have ( [4.36| ) with 
Fi{x,t) = -£^F{x,t). This identity implies ( [4.32| ) with Fi instead of F for 
step functions 

n n 

f = Y.fih-.M^ 9 = Y.9jMc,M^ /i,<?, gC2«. (4.37) 

There is a slight subtlety since T, is not necessarily increasing. However, we 
conclude from ( 4.32| ) and the Parseval equality that for all step functions f,g 



(Fo,Go)i2(K) = (/,5)l2(m+,c2") + / / fixyF{x,t)git)dxdt. (4.38) 

Jo Jo 



INVERSE PROBLEM FOR FIRST ORDER SYSTEMS 



19 



Since a is increasing the assertion now follows from the denseness of the step 
functions in L^Qjjjp(]R+, C^"). To complete the proof it remains to note that the 
equality F{x,t) = Fi{x,t) is a consequence of ( p2|) and (|438|) . □ 



Combining Propositions 4.3 and 4.6 one immediately obtains the following 
theorem. 



Theorem 4.8. Assume that the system (3^) is of class (Tb)- Let a be its 
spectral measure function and cro(A) = ^B^^X. Then with F defined by ( 4.34 ) 
we have the Gelfand-Levitan equation 



F{x,t) + K{x,t) + K{x,s)F{s,t)ds = 0, t < x. 
Jo 



(4.39) 



Remark 4.9. Note that by Proposition 4.6 2. the function F is continuous 
also on the diagonal. In view of ( 4.21| ) the continuity of F at the diagonal 
implies R{x,x) = R{x,x)*. 

Proof. We present a second proof of the Gelfand-Levitan equation based 
on the formula ( [4.34 ) for F, which is similar to |10| and |^^, Chap. 12]. 



For f,ge L2^jjjp(M+,C2") we consider 



I(/,5) := / r 
JR Jo 



dx / dtf{xyY{x,X)da{\)eQ{t,\)*g{t). 



(4.40) 



Substituting (4.5) for Y we find using the Parseval equality and Lemma 4.6 



I(/,9) = (/,9) + 



CO roo 



JO 



f{x)*F{x,t)g{t)dxdt 



+ 



dx / dtf{xy / K{x,s)eQ{s,X)dsda{\)eQ{t,\)*g{t), (4.41) 



ii{f,g)= [ r 

jR Jo 



dx 



ds 



n(/,9) 

f{x)*K{x, s)dxeo{s, X)da{X)eo{t, X)*g{t). 



Writing da = dT, + dao and using Lemma |4.6| we find 

n(/,5) = 



oo roo 

f{x)*K{x,t)g{t)dxdt 
Jo 

oo roo roo 



+ 



Jo 



hence 
I(/,9) = (/,9) + 



oo roo 



Jo 



f{xyK{x, s)F{s, t)g{t)dxdtds, (4.42) 



(4.43) 



f{xy[F{x,y) + K{x,y)+ / K{x,t)F{t,y)dt]g{y)dxdy. 

Jo 

Now if supp/ C [b, oo), suppr? C [0, a], a < 6, then (/,<?) = and 

) 

60(2;, Xy g{x)dx 
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is the c^/f-^Q-transform of 

g{x)+ Rit,x)*g{t)dt 

J X 

which also has support in [0, a], hence by the Parseval equahty l{f,g) = 0. This 
imphes the assertion. □ 

5. The inverse problem 
5.1. The main result. 

Proposition 5.1. Let B = diag(i?i, —B2) he an arbitrary nonsingular self- 
adjoint matrix of signature 0. Let (t(A) be a n x n matrix function satisfying: 

1. If g e Ll^p{R+,C'n and if 

[ Go(A)*dfT(A)Go(A) =0, 

JR 

where Gq is the J^h,o ^transform of g, then g = 0. 

2. The function 



Fix,t):= / eo(x,A)dS(A)eo(t,A)* (5.1) 
with S = o" — (To exists, and has a continuous mixed second derivative 

Then the Gelfand-Levitan equation (4.39) has a unique continuous solution 
K -.n — > M(2n,C). 

Moreover, if F(x,t) is continuously differentiahle, then so is K{x,t). 

Proof. Since for fixed x equation ( 4.39|) is a Fredholm equation it suffices 
to show that the dual equation 

k{t)+ [ k{s)F{t,s)*ds = 0, (5.2) 
Jo 

where k : [0, x] M(2n, C) is square integrable, has only the zero solution. 
Looking at the individual columns in (|5.2D it suffices to show that 

gity + r gisYFit, sYds = 0, ge L\[0, x], C^") (5.3) 



implies g = 0. Extending (7 by to M+ we may consider g as an element of 

^comp ( 



L2 r]R+,C2") and (jsj) implies in view of 2. and Proposition 3. 



= Ibf + / / g{srF{s,t)g{t)dsdt 
Jo Jo 

= \\gf + [ Go(A)*c?E(A)G'o(A) = / Go(A)*da(A)Go(A) 
Jr Jr 

and thus 5 = by 1. 

The proof of C^-smoothness of K{x,t) is similar to that used in [23] and 



10 and is omitted. □ 
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Theorem 5.2. Let B = A\a,g{Bi, —B2) he an arbitrary nonsingular self- 



adjoint matrix of signature as in (1.1). Let 0"(A) be an increasing (right- 
continuous, o"(0) = 0) n X n matrix function satisfying the conditions 1. and 2. 
of Proposition |5.1| . 

Then there exists a unique continuous 2n x 2n matrix potential Q satisfying 
(^) such that the corresponding system (|3.2|) is of class (Tb) and such that a is 
its spectral measure function. Q{x) hasp continuous derivative iff DxD^F{x,t) 
is continuous. 

Conversely, if a is the spectral measure function of the boundary value prob- 
lem (|3.2|) of class (Tb) then the conditions 1. and 2. of Proposition ^]l| hold. 



Proof. The necessity was proved in Lemma 4^ and Proposition i.t . 

To prove the sufficiency we assume that the conditions 1. and 2. of Propo- 
sition |5.1| hold: _ 

i) Starting with ct(A) we define F,F by ( ^.1| ) and ( |5.1| '). Then we consider 
the Gelfand-Levitan equation ( [4.39 ) 



^{x,t) := F{x,t) + K{x,t) + K{x,s)F{s,t)ds = 0, x > t. 

Jo 



(5.4) 



By Proposition 5.1 this equation has a unique continuous solution K : Q ^ 
M(2n,C). 

Then F also equals the right hand side of ( [4.21| ): namely, starting with 
K we consider the operator R of the form ( [4. IT]) and introduce Fi by ( 4.21j ). 
According to Proposition 4.3 Fi and K are connected by equation ( |4.24D. Thus 
F defined by (^) and Fi defined by ( |4.21| ) satisfy the equation ( ^.4| ) and 
therefore we infer from Proposition 4.3 that F = Fi. 

We collect further properties of F: in view of (4.21) we have 



F{x,t) = F{t,x) 



(5.5a) 



By continuity, the equation ( 4.21 ) also holds for x = t and consequently 
R{x,x) is self-adjoint. Therefore, so is K{x,x) = —R{x,x). Furthermore, 



dtF{x,t)B = -BdccF{x,t), 



(5.5b) 



where this equality holds in the distributional sense if F is only continuous. To 
see this let f,g £ C^((0, 00), C^"). In view of ( |43^ ) and (|3lO| ) applied with 
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Q = we calculate 



oo roo 



Jo 



f{xydtF{x,t)Bg{t)dxdt 

f oo roo 



-i I I fix)* F{x,t)-Bdtg{t)dxdt 
Jo Jo ^ 

-i [ {^H,of){\rd^{\)\{^Hfi9)W 

JR 

-i [ {^H,o-Br){\rdT.{\){^H,og){\) 

Jm I 



oo roo 



Jo 



f{x)*Bd^F{x,t)g{t)dxdt. 



Moreover, it follows from ( |5.1| ) and ( 2.11] ) that with some matrix function T[t) 
we have 

mi) = (^)ni)- (5.5c) 

We now define (cf. (p^) 

y(x, A) = eo(x, A) + / K{x,t)eo{t,X)dt (5.6) 



Jo 

and we will show that the properties ( ^.5a[ -c) imply that Y{x, A) satisfies the 
initial value problem 

B-^^^^^ + Q{x)Y{x,\)=\Y{x,X), y(0,A)=eo(0,A)= (5.7) 

where 

Q(x) := iBK{x, x) - iK{x, x)B. (5.8) 

Note that since K{x,x) is self-adjoint Q{x) is self-adjoint, too. Moreover, from 
( ^ ) we also conclude that Q{x) is off-diagonal, i.e. Qa = 0. 



It follows from (|5.5cD that 



F{x, 0)BF{0, t) = T{x)* [Bi - H*B2H]T{t) = 0. (5.9) 



Plugging (|5.9|) into the Gelfand-Levitan equation (|5.4| ) gives 



ir(x,0)5( ^ ) = for xG[0,oo). (5.10) 
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ii) For the moment we assume in addition that F is continuously differen- 
tiable. Then by Proposition 5.1 K also is continuously differentiable. Differen- 
tiating (^.41) we obtain 



Bd^^{x,t) = Bd^F{x,t) + Bd^K{x,t) + BK{x,x)F{x,t) 

/•X 

+ / Bd:^K{x,s)F{s,t)ds = 0, 
Jo 

dt^{x,t)B = dtF{x,t)B + dtK{x,t)B 

/•X 

+ / K{x,s)dtF{s,t)Bds = 0. 
Jo 

Integrating by parts and using ( 5.5b|) and ( 5.10| ) we obtain 



(5.11) 
(5.12) 



K{x,s)dtF{s,t)Bds 



K{x,s)BdsF{s,t)ds 



= / dsK{x,s)BF{s,t)ds - K{x,x)BF{x,t). (5.13) 
Jo 

Adding up ( 5.11 ) and ( ^.12 ) and using ( 5.13| ) and the Gelfand-Levitan 
equation (|5.4|) we obtain 



Bd^K{x, t) + dtK{x, t)B + iQ{x)K{x, t) 

I'X 

+ / {Bd.^K{x,s)^dsK{x,s)B^iQ[x)K{x,s)\F{s,t)ds = ^. 
Jo 

Since the homogeneous integral equation corresponding to the Gelfand-Levitan 
equation ( p.4|) has only the trivial solution (see the proof of Proposition 5.1) we 
infer from (p^) that 



Bd^K{x, t) + dtK{x, t)B + iQ{x)K{x, t) = 0. 



(5.14) 



Since K satisfies the relations ( 5.10 ), ( |5.8D and ( 5.14 ) it follows from Theorem 



4.1 that Y{x, A) (cf. (|5.6| )) satisfies the initial value problem (| 

iii) We now assume that F is just continuous. Assume for the moment that 
for 5 > we have a continuously differentiable matrix function F^ : — > 
M(2n, C) with the properties: 

F^ converges to F as (5 — > uniformly on compact subsets of M^. (5.15a) 
F^ satisfies (|5^ c). (5.15b) 

We fix xq > 0. For < x < xq let Tp be the integral operator in 
C([0,x],C2") defined by (TpDit) := £ f{s)F{s,t)ds. The proof of Propo- 
sition |5.1| shows that —1 spec Tp. Thus for 5 < So{xo) we have —1 spec Tpi 
and the Gelfand-Levitan equations 

{{I + Tps)Ks{x,.)){t) = Ks{x,t) + [ K5{x,s)F^{s,t)ds = -F^{x,t) (5.16) 

Jo 

have (for each fixed x G (0, xq]) unique solutions Ks{x, t), (x, t) G [0, xq]^, which 
converge to as 5 ^ uniformly on [0, xq]^. Since F^ is it can be shown 
(cf. the proof of Proposition 5.1) that Ks is C^, too. 
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Moreover, Ks satisfies ( |5.10| ) for < x < xq which follows fr om ([5.15a ) 
and ( |5.16| ). Now part ii) of this proof shows that Ks also satisfies (|5.14| ) with 
Qsix) := iBKs{x,x) — iKs{x,x)B,x G [0,xo]. Hence, Ks satisfies ( [4.6a| -c) (on 
[0,Xo]^) and therefore, 

/•X 

Ys{x, X) := eo{x, X) + / Ks{x,t)eQ{t, X)dt, < x < xq, 



satisfies the initial value problem { p.7\ ) with Qs instead of Q. 

bmce , x) one concludes as in part i) of this proof that 

QsixT = Qsix). 

Since Ks converges to as (5 — > uniformly on [0,a;o]^, Qs converges to 
Q uniformly on [0, xq]. Thus Y{x,X) satisfies the initial value problem ( ^?^ ) on 
[OjXo]. Since xq was arbitrary Y{x,X) satisfies ( ^.TD on M+. 

It remains to prove the existence of the sequence F^: 

Let F{x,t) := {Fij{x,t))l be the block-matrix representation with re- 
spect to the orthogonal decomposition C^" = 
It follows from and (jsj') that 

Fij{x, t) = fiji^^ix - fijt), fijiO = HigiOH*, (5.17) 

with fii = X^^ , 1 < i < r, and Hi := /„j = /„. Here the map (7 : M — > 
M(n X n, C) is continuous and satisfies g{£,)* = g{—^)- Therefore the maps 
fij : R — > M(nj x nj,C) are continuous and satisfy fijH)* = fji{—C)- We 
note that if the measure S(A) is finite, that is Jjg|(iS(A)| G M(n,C), then 
9{0 = I^e^^^d^iX). 
We put 

A 1 f^^^ 



2<5 J^.s (5.18) 
4(0 := H,g'{i)H* , Ff^{x,t) := /^(/x.x - ^iJt), 

and F'{x,t) := {Ff^{x,t))l^^,. 

Obviously, F^ is continuously differentiable and satisfies (5.15). It is clear 
from KWj that 

. 1 K+s 



Y^j ^9'{s)ds = g\-i), 



(5.19) 



25 

and thus 4(0* = 4(-0- 

In view of ( |5l8| ) and ( p^ F^ satisfies ( ^,b) To prove the property 
( ^ ) for F"^ we note that in view of ([sT?!) and (^.18|) i^^j(0,t) = fij{-fJ.jt) = 
Hig^ {—fijt)H* and consequently 

F^(o,t) = (4(o,t))[,,=i = (if./(-^,t)i/;)i;,.=i =: Qr^(t), (5.20) 

where T^{t) = {g^ {-fiit)Hl g^ {-fi2t)Hl ■ ■ ■ , 9\-f^rt)H;). 
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This proves that satisfies ( |5^ ). Summing up, we have proved that F 
satisfies ( p.l5| ,b). 

iv) Starting with an increasing n x n matrix function cr(A) satisfying the 
conditions 1. and 2. of Proposition we have constructed the boundary value 



problem (|3.2D resp. (5.7). To complete the proof it remains to show that cr{X) 



is, in fact, the spectral function for the problem (5.7). 

Let q{\) be the spectral function of the problem (5/7). Starting with := 



^) — o"o we define Fp by (|5.l|'). Then by Theorem K satisfies the Gelfand- 



Levitan equation ( 4.39| ) with F^. On the other hand, in view of ( |5.4D K satisfies 



the Gelfand-Levitan equation with F instead of Fp. From Proposition 4.2 we 



infer F = Fp. By Remark 4.7 this implies g = a. □ 



Remark 5.3. 1. The case n = 1 and Bi = B2 = 1, i.e. the case of a 
2x2 Dirac system, is due to M. Gasymov and B. Levitan [P], p^, Chap. 12]. 



We note, however, that the proof in |23, Chap. 12] is incomplete, since the 
self-adjointness of Q is not proved. 



2. We also note that following Krein's method [17] L. Sakhnovich |29, 
Chap. 3, §3] has obtained some (implicit) sufficient conditions for a matrix 
measure to be the spectral function of a canonical system. 

5.2. Some complements to the main result. Next we will discuss sev- 



eral other criteria which imply conditions 1. or 2. of Proposition 5.1. For 
brevity, in the sequel we will address them just as "condition 1./2.". As in the 
proof of Theorem ^ we denote by A : ^^(M) ^^(M), (Ag)(A) := \g{\) the 
operator of multiplication by A. Furthermore, we denote by /xj'(Ao) the multi- 
plicity of the spectrum of a self-adjoint operator T at the point Aq- We first 
note some simple facts: 



Remark 5.4. 1. If S(A) is increasing then condition 1. is trivially fulfilled. 
For let g E L2^j^p(]R+, C^") with 

= [Go{\rda{X)Go{X) 

(5.21) 

Go(A)dao(A)Go(A) + / Go(A)(iS(A)Go(A), 



where Go is the ^//^o^transform of g. Since S(A) is assumed to be increasing 
both summands on the right hand side of ( ^.21| ) are nonnegative and hence 0. 
Then Proposition implies g = 0. 

2. Assume that the matrix measure S(A) is finite, i.e. Jjg |(iS(A)| E M(n, C). 
Then condition 2. is obviously fulfilled. 



Recall that a subset X C M is said to have finite density (cf. |21|) if 



limsup4l{2; G ^1 kl < < 00. (5.22) 

Otherwise, X is said to have infinite density. 

Proposition 5.5. Let B = {Bi,—B2) be as in ([Ll|). For an increasing 
n X n matrix function a the condition 
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1'. The set supp^{da) := {A G M | fJ-A^X) = n} has infinite density 
implies condition 1. 

Proof. Let cr(A) = {(7ij{X))'ij^i and g{X) := trcr(A) = crii(A) + ...+(Tnn(A). 
From the inequality 

\(Jij{X) - aij{^)\ < ^ OiiiX) - a,i{iJ.)^Jajj{X) - ajjifi), /i < A, 

we infer that daij{X) is absolutely continuous with respect to dg{X). Hence, by 
the Radon-Nikodym Theorem there exists a density matrix 

$(A) = (<A.,(A))^,=i, (5.23) 

such that 

7(o,A] Ht)dg{t), A > 0, 

Obviously, $(A) > and thus we have 
supp„((i(T) = supp(det $d£>) 

r f^^' 1 (5-24) 

= <^ A G supp{dg) / det ^{X)dg{X) > for all e > ^. 

"/ A — £ 

To check condition 1. let g G L2([0, 6], C^") with 

0=/ Go{Xrda{X)Go{X)= Go{Xy^{X)Go{X)dg{X). 



CT(A) 



Then we have Go (A) = for all A G supp„(d(T). 

On the other hand Go (A) is an entire (vector) function of strict order one 
and hence (cf. ]21[ ) either Go = or the set of its zeros is of finite density. But 
since supp^((ifT) is assumed to have infinite density, we must have G = 0. □ 

Remark 5.6. 1. Condition 1.' of the previous proposition is satisfied if 
supp„((i(T) has at least one finite limit point. 

2. Note that if n = 1 then supp„((i(T) = supp{da) equals the support of the 
Radon measure da. 

Corollary 5.7. Let B he as before and ao{X) = ^B^^X. If 

(i) the measure da is discrete, 

(ii) supp^{da) = supp{dg), 

(iii) the measure dT, = da — daQ, is finite, i.e. Jj^ |(iS(A)| G M(n,C). 
Then a satisfies condition 1.' of Proposition [S^ and hence condition 1. 

Proof. Since d$](A) is finite we have 

a{X) = ^B-^X + C± + o{l), X^±oo, (5.25) 

where G± := S(ibcxD) = G^. Hence 

£»(A) = coA + ci +o(l), A^+oo. (5.26) 
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Since da and hence dg is discrete, we infer from ( |5.26 ) that the set of discon- 
tinuities of Q has infinite density. Hence, by (ii) the set supp„((icT) has infinite 
density. □ 



Remark 5.8. Note that this proof only uses the asymptotic relation (|5.25| ) 

which is slightly weaker than the finiteness of the measure (since the o(l) 
need not be of bounded variation). 

Finally, we give a criterion for the condition 2. 

Proposition 5.9. The condition 2. is fulfilled if the limit 

lim / eo(x,A)dS(A)eo(t,A)* (5.27) 
JlAKA 



exists locally uniformly in x,t. Then indeed F{x,t) is given by ( 5.27 ). 

This is the case if the matrix function 5](A) is integrable with respect to 
Lebesgue measure and satisfies lim S(A) = 0. 

A— >itoo 

Proof. If the limit ( |5.27| ) exists locally uniformly in x,t then we have 
F{x,t) := [ eo(x,A)dS(A)g-o(t,A)* 

JR 

= / / lim / eoix',X)dJ:iX)eo{t',X)*dx'dt' 

Jo Jo ^'^°°J\X\<A 

and we reach the first assertion. 

If S(A) is integrable with respect to Lebesgue measure and satisfies 
lim S(A) = 0, then we apply integration by parts for Lebesgue-Stieltjes in- 

A— >±oo 

tegrals to obtain 

eo(x,A)(iS(A)eo(t, A)* = eo{x, X)T.{X)eo{t, X) I^Z^a 

A|<A 

- / {dxeo{x,X))^{X)eo{t,X)+eo{x,X)^{X){dxeo{t,X))dX. 

J\X\<A 

Since eo(x,A) is uniformly bounded and dxeo{x,X) is uniformly bounded for 
|x| < -R for each R, we reach the conclusion. □ 

Corollary 5.10. The spectral function can be prescribed on an arbitrary 
finite interval. More precisely, there exists a boundary value problem (3^) with 
continuous Q satisfying (O) and such that its spectral function a{X) coincides 
on an arbitrary finite interval with a prescribed increasing nxn spectral measure. 

Proof. This follows from the fact that if S(A) is constant outside a com- 
pact interval then it satisfies condition 2. by Remar k |5.4| (or the previous 
proposition) and it satisfies condition 1. by Proposition |5.5|. □ 



Conjecture 5.11. We conjecture that condition 1. in Theorem is ob- 
solete in general. 
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6. Some generalizations, comments, examples 

6.1. Generalization of the main result Theorem |5.2| . Before we have 
investigated an operator L of the form ( |3.2| ) starting with the operator Lq (with 
(5 = 0). This has an obvious generalization. Namely, we may investigate two 
operators Li := Li^h,I-'2 := ^2,^ and consider L2 as a perturbation of Li. 
More precisely, let 

L,-\b±+Q,. (6.1) 

and 

&{Lj) = {/ e ^{L*) I /2(0) = /?/i(0)}, Bi = H*B2H. (6.2) 

Furthermore, let Yj be the 2n x n matrix solution of the initial value problem 
( ^ ) (with Lj instead of L). If both operators Lj are of class (Tg) then Yj 
admits the representation Yj(.,X) = (I + Kj)eo{., X), where Kj is a Volterra 
operator with kernel Kj{x,t). Therefore 



Y2ix,X) = {iI + K)Yi{.,X))ix)=Yi{x,X)+ Kix,t)Yi{t, X)dt, (6.3) 

Jo 

where 

/ + K = (l + ir2)(/ + -?^i)"^ (6.4) 

Repeating the arguments used in the proof of Theorem one concludes 
that if Qi,Q2 G C\R+,M{2n,C)) then K G (7^(0, M(2n, C)) and, moreover, 
satisfies the following Goursat problem 

BdxK{x, t) + dtK{x, t)B + iQ2ix)K{x, t) - iK{x, t)Qi{t) = 0, (6.5a) 
BK{x, x) - K{x, x)B = i{Qi{x) - ^2(2;)), (6.5b) 

K{x,Q)b(^^=Q. (6.5c) 

We also note that ( |6.5a| )- (|6.5c| ) may be deduced directly from (6.4) and ( [4.6| ) 
for Ki,K2. For example ( |6.5bD follows from ( [4.6b ) and the identity K{x,x) = 
K2{x,x) - Ki{x,x). 

Putting R:= {I + K)"^ - I we obtain from (O) 



Yi{x,X) = {{I + R)Y2{.,X)){x)=Y2{x,X)+ R{x,t)Y2{t, X)dt. (6.6) 

Jo 

Since Proposition |4.3| remains valid in the case under consideration, the 



following result, being a complete analog of Proposition |4.6| , may be obtained 
in the same way as Proposition 

Proposition 6.1. Let crj{X) be the nxn spectral function (cf. Proposition 
3^) of the operator Lj,j = 1,2, and S := (J2 — ai. 

1. Let Lj be of class (Tb) and let F{x,t) be defined by ( 4.31| ) with R{x,t) 
being the kernel of the transformation operator (|6.6|). Then we have for all 
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Fi(A)*dS(A)Gi(A) = / / f{xyF{x,t)git)dxdt, (6.7) 
Jo Jo 

where Fi and Gi are the ^h,Qi -transforms of f and g respectively. 
2. Again assuming Lj to be of class (Tb) we put 



Yi{x,X) := / Yi{t,X)dt. 
Jo 

Then the function 

F{x,t):= [ Yi{x,X)di:{X)Yi{t,X)* (6.8) 



exists and has a_ continuous mixed second derivative which coincides with 
Fix,t), i.e. ^,Fix,t) = F{x,t). 

3. Conversely, given any increasing n x n matrix function 02 put S := 
(72 — (Ji. If the integral (|6.8| ) exists and has a continuous mixed second derivative 
Fi{x,t) := ^F(x,t) then holds for all f,g e L2^j„p(M+, C^") with Fi 

instead of F. 

Again, we emphasize that 3. holds for arbitrary L of the form (|3.2| ) not 
necessarily being of class {Tb). 



Combining Propositions |6.1| and |4.3| we arrive at the Gelfand-Levitan equa- 
tion: 

Proposition 6.2. LetLj be of class (Tb) and letaj be the spectral function 
of the problem ( |3.2D with Qj = j = 1, 2, instead of Q. Then with F defined 
by ( |6.8D we have the Gelfand-Levitan equation ( [4.39] ). 

Now we are ready to present a generalization of the main result (Theorem 



Theorem 6.3. Let o"i(A) be the spectral function of the operator Li of the 
form ( [6.11 ). For an increasing n x n matrix function o"(A) to be the spectral 
function of the boundary value problem ( p.2[ ) with (unique) continuous 2n x 2n 
matrix potential Q satisfying ( |4.3| ) it is sufficient that the following conditions 
hold: 

1. //5Gi?omp(K+,C2") and if 

Gi{X)*da{X)Gi{X) = 0, 

where Gi is the ^h,Qi ^transform of g, then g = 0. 

2. The function 

F{x,t) := [ Yi{x,X)d^{X)Yi{t,X)*, Yi{x,X):=[ Yi{t,X)dt, (6.9) 
Jr Jo 

with S = cr — o"! exists and has a continuous mixed second derivative 



30 



MATTHIAS LESCH AND MARK MALAMUD 



Moreover Q has m continuous derivatives if and only if D^DY^F{x, t) exists 
and is continuous. 

Again, if we content ourselves to operators of class (Tb) then the conditions 
1. and 2. are also necessary. 

Sketch of Proof. The necessity is proved in just the same way as Lemma 



|4.5| and Proposition 4.6 



Sufficiency: Starting with (t(A) we define F,F by ( |6.8D with S(A) := a{X) 
cri(A). Then we consider the Gelfand-Levitan equation 



F{x,t) + K{x,t) + K{x,s)F{s,t)ds = 0, t < x. 
Jo 



(6.10) 



with F defined by (|6.§| ). Following the proof of Proposition 5J one concludes 
that ( 6.1C| ) has a continuous solution K : Q ^ M(2n, C). Next we define Y{x, A) 
setting Y{., A) = {I + K)Yi{., A) and show that Y{x, A) satisfies the initial value 
problem (|5.7D with 

Q{x) = Qi{x) + iBK{x, x) - iK{x, x)B. (6.11) 

Since Qi satisfies (|4.3| ) we infer from ( |6.11 ) that Q also satisfies ([4.3|). Moreover 
the self-adjointness of Q may be proved as in the proof of Theorem |5.1| . 

Furthermore, we note that if F is continuously differentiable it satisfies the 
equality 

BD^F{x,t) + DtF{x,t)B = -iQi{x)F{x,t) + iF{x,t)Qi{t) (6.12) 



and according to Proposition |5.1| K is continuously differentiable, too. If F is 
just continuous then (6.12) still holds in the distributional sense. This is shown 
similar to ( |5.5bD . 

Since Yi(0,A) = (^) we may argue exactly as in ( 5.5cj ), ( |5.9[ ), ( ^.10 ) to 
obtain 

K{x,0)b(^^^ =0, for xe[0,oo). (6.13) 

In view of ( |6.10D -( p.l3| ) the relation (6.5a) for K is proved along the same 
lines as part ii) of the proof of Theorem ^.2|. 

Thus K satisfies the initial value problem ( |6.5a| )- (l6.5c ). Therefore Y{x,X) 
satisfies the initial value problem ( ^.7]) with Q defined by ( |6.1lD . 

If now F is just continuous then one proceeds as in part iii) of the proof of 
Theorem 5.2. 



That a is indeed the spectral function of the problem (5/7) with Q from 
( 5.11 ) is shown as part iv) of Theorem 5^2. Instead of (|5.1|), T heorem |4^ , ( ^.4| ) 
one uses (^ 



and Proposition 



Proposition 5.2, (p.l 



6.1 



and Proposition 
□ 



6.2. The degenerate Gelfand Levitan equation. We discuss solutions 
of the Gelfand-Levitan equation in the special case where 5](A) is a step func- 
tion: 

We consider the situation of Theorem 6.3 and fix an operator Li of the form 
(|0]) with spectral function cri(A). 
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Let A G M(n, C) be a hermitian nonnegative matrix and 

S(A) :=^l[„,oo)(A) (6.14) 

an increasing step function with one jump of "height" A. 
We show that 

0- := 0-1 + E (6.15) 
is the spectral function of the boundary value problem (p. 21) for some (unique) 



continuous self-adjoint 2n x 2n-matrix potential Q satisfying ( |4.3| ). 

Since jumps of the spectral function correspond to eigenvalues this shows 
in particular that for a given potential Qi and given real number a there is a 
potential Q such that 

spec(Li + Q — Qi) = spec(Li) U {a}. (6.16) 

For the proof we have to verify the conditions 1. and 2. of Theorem 5.3. 
By Remark 5.4 condition 1. is fulfilled since A is nonnegative. To verify 2. we 
calculate 



F{x,t) = / Yi{x,X)dJ:{X)Yi{t,X)* 
Jr 

= Yi{x,a)AYi{t,ay. 
Obviously, this has a continuous mixed second derivative, namely 

F{x,t) := ^^F{x,t) = Yi{x,a)AYi{t,ay. (6.17) 

In this case we can solve the Gelfand-Levitan equation explicitly. First we 
introduce for a; > 

T{x) := [ Yi{s,ayYi{s,a)ds. (6.18) 
Jo 

From 

Yi{0,ayYi{0,a) = I + H*H > I 

we infer that T{x) > is positive definite for x > 0. 
We put for t < X 

K{x,t) := -Yi{x,a)AYi{t,ay + 

+Yi{x, a)AT{x)A^/'^{I + A^/^T{x)A^/'^)-^A^/'^Yi{t, ay 

= -Yi{x,a)A^/'^{I + A^/^T{x)A^/^)-^A^/'^Yi{t,ay. 

(6.19) 

Note that (/ + A^^'^T{x)A'^^'^) > I is positive definite, thus invertible. We 
abbreviate S{x) := A^/'^{I + A^/'^T{x)A^/'^y^A^/'^. U A is positive definite then 
we simply have S{x) = {A~^ + T{x))~^ . 

One immediately checks that K{x,t) solves the Gelfand-Levitan equation 
( 16.101 ) corresponding to F and consequently determines Q by means of of ( 6.11| ). 

Summarizing the previous considerations we arrive at the following propo- 
sition. 
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Proposition 6.4. Let Li be an operator of the form (p^), (|OD with the 



spectral function (Ti(A) and let S(A) be of the form ( |6.14 ). Then a = ai + is 
the spectral function of the boundary value problem (|3.2) with 2n x 2n matrix 
potential 

Q{x) = Qi{x) + i{Yi{x,a)S{x)Y{{x,a)B - BYi{x,a)S{x)Y{{x,a)}. 



Corollary 6.5. Under the assumptions of the previous Proposition p.4\ let 
Qi = (i.e. Li = —iB-^). Then the 2n x 2n matrix potential corresponding 
to the spectral function cr(A) = 2^x^^In + ^^(A) with one jump of "height" A is 
given by 

Q{x) = ie'''^''^S{x)B - BS{x)}e~'^^'''', 
where S{x) := {lj)A^I'^{I + xA^/^{I + H*H)A^/^)-^A^I'^{I,H*). 

If E is a general increasing step function then the Gelfand-Levitan equation 
is still solvable. However, we do not have such an explicit formula: 

Proposition 6.6. Let Li be an operator of the form ( |6.1| ), ( |6.2D with spec- 
tral function (Ti(A). Furthermore, let -co < ai < . . . < Ur < oo be real numbers 
and Aj £ M(n, C) nonnegative matrices. 

Then for the increasing step function 

r 



there exists a unique continuous matrix potential Q satisfying (O) such that 
ai + T, is the spectral function of Li + Q — Qi. 



Namely, Q is uniquely determined by (6.11) with K{x,t) being the solution 



of the Gelfand-Levitan equation (6.10) 



Proof. This follows by induction from Theorem 6^ and the preceding 
discussion. 



The conditions 1. and 2. of Theorem 6.3 can also immediately be checked 
directly: namely, condition 1. is fulfilled in view of Remark |5.4| since S is 
increasing. Condition 2. follows immediately from 

r r 

F{x,t) ='^Yi{x,aj)AjYi{t,aj)* and F{x,t) ='^Yi{x,( 2j)AjYi{t, Oj)*. 

□ 

6.3. On unitary invariants of 2n x 2n systems. It is well-known that a 
selfadjoint operator A in a Hilbert space is uniquely determined (up to unitary 
equivalence) by the spectral type [Ea] and the multiplicity function Nej^. In this 
section we will show that there exist potentials Q such that the corresponding 
operator Lh has constant multiplicity one and [E] is of pure type (absolute 
continuous, singular continuous, pure point). 

Definition 6.7. An increasing function /x : M ^ M on the real line will be 
called p-admissible if there exists a strictly increasing sequence of real numbers, 
{xu)u€Z, such that 
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1. Xq = 0, 

2. the sequence {x^+i — Xy) is square summable, 

3. lim Xy = zbco, 

4. n{ayj) < /i(ai/j+i), where a^j := x^ + j2^'"-^P{xi^+i - x^), 
i = 0,... ,2"+f-l. 

In particular, a strictly increasing function /i is p-admissible for any p. 
We will show that for a p-admissible increasing function /i there exists an 
operator Lu of the form ( |3.2D such that its spectral measure E := E^^ satisfies 

[E] = [d/u], Ne{x) = p for /i- a.e. x G M. 

In particular there exist 2n x 2n systems such that each point in the spec- 
trum has multiplicity one. To prove this result we will use the criteria from the 
end of Section 5. 

Proposition 6.8. Let B = {Bi,-B2) G M(2n,C) be a matrix as in (|]|) 
and H G M(?7-,C) as in ( p. 8] ). Let fi be a p-admissible increasing function on 
the real line, 1 < p < n. Then there exists a continuous potential Q satisfying 
( [4.3| ) and such that the corresponding operator Lh is unitary equivalent to the 
operator Ap = ©fAi , where 

Ai : l2(M) ^ LliR), Ai/(A) = A/(A). 

Proof. Let (ipj) be the Rademacher functions Sec. 1.3], i.e. ipi ■.'R ^ 
M is a function of period one, such that 

Mx) = l I' i^'^f (6.20) 
[-1, i < X < 1, 

and 

iPj{x)=M'^'^)- (6.21) 



(In [31, Sec. 1.3] one puts tjji{l/2) = 0). ipj takes values ±1. The set (ipj) is 
orthonormal in L^[0, 1] and 

V'j(x)dx = 0. (6.22) 



We put 

V'o(A) 



/V'i(A) ^2(A) ... V'n(A) 

M^) V3(A) ... Vn+l(A) 
\ V'p(A) '0p+l(A) ... Ipn+p-lW J 



(6.23) 



V^(A) := M ^ ), for X, < A < x,+i, (6.24) 

^(A) := p-ii^Sr'/'V'(A)*V(A)i?r'^'- (6-25) 
<I>(A) is a symmetric nonnegative matrix of rank p for each A G M. 
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Furthermore, let 



fix) := — f — -, auj <x< au,j+i, (6.26) 



and put 



(X) =U\o,x]f(t)Mt), A>0, 
\-Iix,o]f(t)dKt), A<0. 



(6.27) 



(6.31) 



Obviously the measures dg and dfj, are mutually equivalent and 

Q{a^j) = a^j- (6.28) 

Finally we put 

f L ,1 ^(t)dg(t), A > 0, 

a(X) := <^ ^^V^ . (6.29) 

l-/(A,o]^W^^'(*)' ^<0- 
Note that in view of ( |6.22D and the orthonormality of the Rademacher 
functions we have cr{x^) = ^B^^Xu. Again, by the orthonormality of the 
Rademacher functions and the fact that the entries of ip are constant on the 
intervals (oi,j, Oi/j+i], we have for x^ < X < x^j^i 

\\cj{X) - ^B^^x.W < [ \mt)\\dQit) < Cix,+i - X,), (6.30) 

thus 

/ ll^(A) - ^B^^HdX <C' (x^+i - x^)dX 

= C'{x^+i - Xyf. 
Since {x^^i — x^)^ is square summable, we infer that the function 

S(A) :=a(A)-^7o(A) 
is integrable with respect to Lebesgue measure and lim S(A) = 0. In view of 

A— >itoo 

Proposition |5.9| it satisfies condition 2. 

To show that it satisfies condition 1. let g G L^([0, 6], C^") with 

/OO /'OO 
Go(A)*d^7(A)Go(A) = / G^{Xr^{X)Go{X)dQ{X). 
-OO J — OO 

In view of 4. of the definition of admissibility we infer that in each inter- 
val [xy^Xy+i] there exist points Xyj,j = 0, ..,2"+?' — 1 , such that the vector 

—1/2 

B^ Go{Xuj) lies in the null space of the matrix ijjQ{j2~^~P) . Since (x,y+i —x,^) 
is square summable, each of the sequences {Xuj)u is a sequence of infinite den- 
sity. Noting that Go (A) is an entire (vector) function of strict order one and of 

I I —1/2 

finite type we infer as in the proof of Proposition 5.5 that B^ Go (A) lies in 
the null space of the matrix ipQ(j2~^~P), j = 0, ...,2"'~^'P — 1, for each A. It is 
easy to check that the intersection of these null spaces is 0, hence Go (A) = 0. 



By Theorem 5.2 there exists Q satisfying ( [4. 3D such that a is the spectral 
measure function of the corresponding operator Lh- This proves the theorem. 
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□ 

Corollary 6.9. For 1 < p < n there exist continuous potentials Q : 
M+ — > M(n, C) satisfying ( |4.3| ) such that the corresponding operator Lh has 

i) absolute continuous spectrum of multiplicity p, 

ii) singular continuous spectrum of multiplicity p, 

iii) pure point spectrum of multiplicity p. 

Proof. We only have to note that there exist admissible increasing func- 
tions /^ac, A'sc) A^pp such that the measures d/^ac, dA*sc, c?/^pp are absolute contin- 
uous, singular continuous, discrete, respectively. Such measures obviously ex- 
ist. □ 
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